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We compute the pion inclusive momentum distribution in a heavy-ion collision, assuming ther- 
mal equilibrium and accounting for boundary effects at the time of decoupling. We calculate the 
chemical potential corresponding to an average pion multiplicity in central collisions and explore the 
consequences of having the pion system produced close to the critical temperature for Bose-Einstein 
condensation. 
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*y. , In recent years, much experimental effort has been devoted to the production of matter at high densities and 

y—( ' temperatures in relativistic nuclear collisions. A main goal of such experiments is to detect the transition from 

_ , hadronic matter to the quark gluon plasma (QGP). The expectation is that by colliding heavy systems, as opposed 

to single hadron-hadron collisions, the chances of producing a highly compressed state of matter, such as the QGP, 

are increased. Therefore, any experimental signal for the production of QGP is better distinguished if we compare 

^^ ■ similar signatures to those obtained from the collisions of smaller systems such as nucleon-nucleon (n-n) collisions. 

In this spirit, recent experiments have concentrated their efforts on measuring one of those signatures, namely, the 
I . inclusive single particle momentum distributions for some of the most abundantly produced particles after a heavy 
£> ion reaction such as pions. 

The most striking feature reported [pj — fa] is an enhancement at low as well as at high transverse momentum on 
the inclusive single pion distributions as compared to n-n collisions. This observation has sparked a great deal of 
theoretical effort trying to explain the origin of such peculiar behavior [pi- 113]. Li and Bauer M have supported the 
idea that a superposition of two Maxwell-Bolzmann (M-B) distributions with different temperatures might account 
for the observed behavior. Atwater et al M and Lee and Heinz M have suggested the possibility of transverse 
flow. Kataja and Ruuskanen [H]] have fitted a non-zero chemical potential to the Bose-Einstein (B-E) distribution. 
Mostafa and Wong [pT| stressed the importance of properly accounting for boundary effects at freezout. Sollfrank et 
al [1J] and G.E. Brown el al |13| (see also ref. g) included the influence of resonant decays into the picture. The 
effects of Coulomb final state interactions in describing the momentum distribution of charged particles has also been 
Ph! pointed out (T3J. 

While the expansion, boundary and Coulomb effects have been amply discussed in the literature, one of the main 
features of systems described by B-E statistics has not been given enough consideration, namely, the fact that close 
to some temperature T c , bosons occupy predominantly the lowest energy state [p.5| , flq| . This tendency of bosons to 
bunch together leads to B-E condensation at T c and could be responsible for the enhancement of the distribution at 
low pt- It is remarkable, as we will latter show, that the pion multiplicities reached in some heavy ion experiments, 
together with the size of the pion source, as measured for example by Hanbury-Brown Twiss (HBT) interferometry, 
yield a value for T c close to the temperatures which are consistent with typical abundances fllTf . 

On the other hand, when the system of pions can be treated as being confined just before freezout and the wave 
functions for the states satisfy a given condition at the confining boundary |ll| , the states form a discrete set. In this 
case, the density of states contributing to the momentum distribution is larger at high p t compared to a simple B-E 
distribution. This effect, together with an initial transverse flow at freezout, could account for the enhancement of 
the distribution at high p t . 

In this paper, we compute the momentum distribution for pions produced in a relativistic heavy ion collision, 
assuming that at freezout, the pions are in thermal equilibrium and incorporating boundary effects. We reserve the 
discussion of possible expansion effects for a following up work. 

Recall that a system which consists of a fixed number N of weakly interacting, spin zero bosons, in thermal 
equilibrium at a temperature T is described by a grand-canonical ensemble obeying B-E statistics. If Ei represents 
the energy of a single particle state, labeled by i, then the number of particles, chemical potential p, and temperature 
are related by 



o 

:a 

Oh 

6. 

> 

X 



'Current address: SISSA, Via Beirut 4, 340f4, Trieste Italy 



N = E e i E J )/T _ i • (i) 

i 

Since n(Ei) cannot be negative, we have the condition that /i < E , where E is the lowest single particle energy 
state. For a fixed number of particles, p is a function of T. The peculiar characteristic of eq. (fil) is that when T 
approaches a certain critical value T Cl /i approaches the value E and thus, this energy state becomes more populated 
as T gets closer to T c . 

One can estimate the value of T c by considering the continuum limit of eq. (|l|). The result is that for a system of 
weakly interacting relativistic bosons, T c is given implicitly by 
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where V is the volume, m is the boson's mass and K 2 is the modified Bessel function of the second kind and order 2. 
However, the above limit assigns a weight zero to the lowest energy state. Then it is clear that when T — > T c , one needs 
to return to the discrete picture and account for the contribution of the individual energy states to quantities such as 
the momentum distribution. For this purpose, we are required to make general assumptions about the evolution of 
the system. 

We assume that the system of pions (of a particular species) produced after a heavy ion collision is in thermal 
equilibrium at its time of formation. This corresponds to assuming that the total collision rate before the time of 
decoupling is high compared, for instance, to the expansion rate. 

For a given average multiplicity, we can then use a grand-canonical ensemble to describe the statistical properties 
of the pion system. But statistics alone is not enough to describe the system's evolution, we also have to account 
for the fact that at the time of decoupling, the pion system has a finite size and is confined within the boundaries 
of a given volume. The shape of this volume is certainly an important issue. Bjorken dynamics Eg ], for instance, 
implies an overall longitudinal expansion and thus that cylindrical geometry is better suited. Moreover, the time of 
decoupling is not necessarily the same over the entire volume J19[ . However, in order to gain physical insight into the 
problem, we will study the case in which the confining volume has spherical shape and the decoupling time is unique 
in the cm. frame. We thus consider the following scenario: At the time of decoupling, when strong interactions have 
ceased, the system of pions (of a given species) is in thermal equilibrium and is confined within a sphere of radius R 
(fireball) as viewed from the center of mass of the colliding nuclei. 

We start by estimating the critical temperature for the onset of B-E condensation. Fig.[j] shows plots of N versus 
T c computed from eq. (0) for different values of R. Recall that the total multiplicity in a heavy ion reaction is a 
function of the invariant energy y/s in the collision. The multiplicity increases logarithmically with y/s. At AGS 
energies (y/s ~ 5 A GeV) for example, the average pion multiplicity per event produced in central collisions is on the 
order of 400-500 and thus the number of pions of a particular kind is roughly a third of the above. From Figjlj, we 
notice that the value of T c for a number of pions of a particular species between 100-200 is fairly high and decreases 
when the volume increases. From the discussion above, the contribution from the discrete energy states has to be 
properly accounted for. Let us proceed to compute this contribution. We first solve for the relativistic wave function 
corresponding to the stationary states of a free particle inside a sphere of radius R 

^2 -V 2 +m 2 ^(r,i) = 0. (3) 

We impose the boundary conditions corresponding to a rigid sphere, namely 

i>(R,t)=0, (4) 

to describe the initial particle confinement at freezout. The normalized solutions to eq. (0), together with the boundary 
condition, eq. (Q), are |20) 
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where J„ is a Bessel function of the first kind and lj m / (f) are the spherical harmonics. The quantum number k is 
given by the solution to 



Ji+i/2(kR) = . (6) 



The energy eigenvalues are related to k by 



E kl = ^/k 2 + to 2 . (7) 



The normalized contribution to the momentum distribution from the energy state with quantum numbers k,l,m' is 
proportional to the absolute value squared of the space Fourier transform of eq. (o) 

4>klm'if) = (22?fc|)|^klm' (f)\ 2 , (8) 

where 

V*w (# = t^W / ^re-V^idm' (r) ■ (9) 

Due to the azimuthal symmetry of the problem, the wave function in momentum space does not depend on the 
quantum number m' and is a function only of the momentum magnitude 

Ipklm' if) = Ipkl (p)S m '0 ■ (10) 

Accordingly, the thermal momentum distribution is given by 
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where <f>ki (p) is given explicitly by 

<t>ki(p) 
and (i is computed from 
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for fixed N. In general, we can expect an enhancement at high transverse momentum, relative to a simple B-E 
distribution, due to the finite size of the system at freezout, which results in a higher density of states contributing 
at high p than at low p. The shape of the distribution at low p should also be affected by a finite chemical potential. 

Fig.@ shows the transverse momentum distribution md 2 N / m 2 dmtdy at fixed y(= 1.4), as a function of pt calculated 
from eqs. (|ll| ) and (^2|) for several values of R, T and N and for /i computed from eq. (|l3|). For the chosen set of 
parameters, the main effect is the bending upwards of the distribution at high p t relative to a simple exponential with 
the same temperature. We also find, in agreement with reference [|ll| , that the distribution starts deviating from a 
simple exponential fall off at smaller intermediate values of p t as R decreases. 

Fig.y shows the distribution for R = 8 fm, T = 150 MeV, N — 150, compared to a simple B-E distribution with the 
same parameters and a chemical potential corresponding to the same number of particles. Notice how the distribution 
in terms of discrete states deviates from the simple B-E distribution at high momentum and that, since the parameters 
are far from the critical region for B-E condensation, both distributions coincide low pt- The situation changes at 
low p when the parametres are close to the critical region, this is depicted in figJU where we show the distributions 
for R = 6 fm, T = 120 MeV, TV = 200. The difference stems from the way the condensate contribution is included 
in the discrete states picture compared to the continuos case. In the former, the condensate contribution is spread 
over different states, although it is mainly concentrated in the lowest ones; in the latter, the condensate contribution 
is restricted to the lowest energy state, p = 0. 

In summary, we propose to compute the pion inclusive distribution accounting for density effects at decoupling by 
calculating the chemical potential corresponding to a given multiplicity. Possible boundary effects can be included by 
a description in terms of a set of discrete states. If the pion system is produced close to the critical region for B-E 
condensation, a description in terms of a discrete set of states lends itself for the inclusion of the contribution from 
the lowest energy states. We also emphasize that the condensation is a high density phenomenon and thus depends o 
both, the pion multiplicity and the freezout volume. The presence or absence of a condensate could be used to indicate 
for instance the degree of transparency in a central collision and have consequences for correlation experiments. 
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Figure Captions 
FIG. 1. Critical temperature vs. Number of particles for different values of the radius R at decoupling. 

FIG. 2. Invariant momentum distribution rnd 2 N/m 2 dm t dy at y — 1.4, for several values of the parameters R, T, N and 
with the computed values of /j,. Notice the bending upwards of the distribution at high p t and the shape at low p t . 

FIG. 3. Invariant momentum distribution md 2 N / m 2 dm t dy at y — 1.4, for R = 8fm, T = 150MeV and N = 150 corresponding 
to a value of /i = 94.3MeV. Shown also is a simple B-E distribution with the same set of parameters corresponding to a value 
of n = 67.4MeV. 

FIG. 4. Invariant momentum distribution md 2 N / rn 2 dm t dy at y — 1.4, for R = 6fm, T = 120MeV and N = 200 corresponding 
to a value of /i = 169.1MeV. Shown also is a simple B-E distribution with the same set of parameters corresponding to a value 
of n = 83.5MeV. 
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